Search for: All records

Galactic nuclei showing recurrent phases of activity and quiescence have recently been discovered. Some have recurrence times as short as a few hours to a day and are known as quasi-periodic X-ray eruption (QPE) sources. Others have recurrence times as long as hundreds to a thousand days and are called repeating nuclear transients. Here we present a multiwavelength overview of Swift J023017.0+283603 (hereafter Swift J0230+28), a source from which repeating and quasi-periodic X-ray flares are emitted from the nucleus of a previously unremarkable galaxy at ∼165 Mpc. It has a recurrence time of approximately 22 days, an intermediary timescale between known repeating nuclear transients and QPE sources. The source also shows transient radio emission, likely associated with the X-ray emission. Such recurrent soft X-ray eruptions, with no accompanying ultraviolet or optical emission, are strikingly similar to QPE sources. However, in addition to having a recurrence time that is ∼25 times longer than the longest-known QPE source, Swift J0230+28’s eruptions exhibit somewhat distinct shapes and temperature evolution compared to the known QPE sources. Scenarios involving extreme mass ratio inspirals are favoured over disk instability models. The source reveals an unexplored timescale for repeating extragalactic transients and highlights the need for a wide-field, time-domain X-ray mission to explore the parameter space of recurring X-ray transients. 

Abstract Pulsar timing array experiments have recently uncovered evidence for a nanohertz gravitational wave background by precisely timing an ensemble of millisecond pulsars. The next significant milestones for these experiments include characterizing the detected background with greater precision, identifying its source(s), and detecting continuous gravitational waves from individual supermassive black hole binaries. To achieve these objectives, generating accurate and precise times of arrival of pulses from pulsar observations is crucial. Incorrect polarization calibration of the observed pulsar profiles may introduce errors in the measured times of arrival. Further, previous studies have demonstrated that robust polarization calibration of pulsar profiles can reduce noise in the pulsar timing data and improve timing solutions. In this paper, we investigate and compare the impact of different polarization calibration methods on pulsar timing precision using three distinct calibration techniques: the Ideal Feed Assumption (IFA), Measurement Equation Modeling (MEM), and Measurement Equation Template Matching (METM). Three NANOGrav pulsars—PSRs J1643−1224, J1744−1134, and J1909−3744—observed with the 800 MHz and 1.5 GHz receivers at the Green Bank Telescope (GBT) are utilized for our analysis. Our findings reveal that all three calibration methods enhance timing precision compared to scenarios where no polarization calibration is performed. Additionally, among the three calibration methods, the IFA approach generally provides the best results for timing analysis of pulsars observed with the GBT receiver system. We attribute the comparatively poorer performance of the MEM and METM methods to potential instabilities in the reference noise diode coupled to the receiver and temporal variations in the profile of the reference pulsar, respectively. 

Abstract Based on the rate of change of its orbital period, PSR J2043+1711 has a substantial peculiar acceleration of 3.5 ± 0.8 mm s^–1yr^–1, which deviates from the acceleration predicted by equilibrium Milky Way (MW) models at a 4σlevel. The magnitude of the peculiar acceleration is too large to be explained by disequilibrium effects of the MW interacting with orbiting dwarf galaxies (∼1 mm s^–1yr^–1), and too small to be caused by period variations due to the pulsar being a redback. We identify and examine two plausible causes for the anomalous acceleration: a stellar flyby, and a long-period orbital companion. We identify a main-sequence star in Gaia DR3 and Pan-STARRS DR2 with the correct mass, distance, and on-sky position to potentially explain the observed peculiar acceleration. However, the star and the pulsar system have substantially different proper motions, indicating that they are not gravitationally bound. However, it is possible that this is an unrelated star that just happens to be located near J2043+1711 along our line of sight (chance probability of 1.6%). Therefore, we also constrain possible orbital parameters for a circumbinary companion in a hierarchical triple system with J2043+1711; the changes in the spindown rate of the pulsar are consistent with an outer object that has an orbital period of 60 kyr, a companion mass of 0.3M_⊙(indicative of a white dwarf or low-mass star), and a semimajor axis of 1900 au. Continued timing and/or future faint optical observations of J2043+1711 may eventually allow us to differentiate between these scenarios. 

Abstract Noise characterization for pulsar-timing applications accounts for interstellar dispersion by assuming a known frequency dependence of the delay it introduces in the times of arrival (TOAs). However, calculations of this delay suffer from misestimations due to other chromatic effects in the observations. The precision in modeling dispersion is dependent on the observed bandwidth. In this work, we calculate the offsets in infinite-frequency TOAs due to misestimations in the modeling of dispersion when using varying bandwidths at the Green Bank Telescope. We use a set of broadband observations of PSR J1643−1224, a pulsar with unusual chromatic timing behavior. We artificially restricted these observations to a narrowband frequency range, then used both the broad- and narrowband data sets to calculate residuals with a timing model that does not account for time variations in the dispersion. By fitting the resulting residuals to a dispersion model and comparing the fits, we quantify the error introduced in the timing parameters due to using a reduced frequency range. Moreover, by calculating the autocovariance function of the parameters, we obtained a characteristic timescale over which the dispersion misestimates are correlated. For PSR J1643−1224, which has one of the highest dispersion measures (DM) in the NANOGrav pulsar timing array, we find that the infinite-frequency TOAs suffer from a systematic offset of ∼22μs due to incomplete frequency sampling, with correlations over about one month. For lower-DM pulsars, the offset is ∼7μs. This error quantification can be used to provide more robust noise modeling in the NANOGrav data, thereby increasing the sensitivity and improving the parameter estimation in gravitational wave searches. 

Abstract We test the impact of an evolving supermassive black hole mass scaling relation (M_BH–M_bulge) on the predictions for the gravitational-wave background (GWB). The observed GWB amplitude is 2–3 times higher than predicted by astrophysically informed models, which suggests the need to revise the assumptions in those models. We compare a semi-analytic model’s ability to reproduce the observed GWB spectrum with a static versus evolving-amplitudeM_BH–M_bulgerelation. We additionally consider the influence of the choice of galaxy stellar mass function (GSMF) on the modeled GWB spectra. Our models are able to reproduce the GWB amplitude with either a large number density of massive galaxies or a positively evolvingM_BH–M_bulgeamplitude (i.e., theM_BH/M_bulgeratio was higher in the past). If we assume that theM_BH–M_bulgeamplitude does not evolve, our models require a GSMF that implies an undetected population of massive galaxies (M_⋆≥ 10¹¹M_⊙atz> 1). When theM_BH–M_bulgeamplitude is allowed to evolve, we can model the GWB spectrum with all fiducial values and anM_BH–M_bulgeamplitude that evolves asα(z) =α₀(1 +z)^1.04±0.5. 

Abstract Pulsar timing array observations have found evidence for an isotropic gravitational-wave background with the Hellings–Downs angular correlations between pulsar pairs. This interpretation hinges on the measured shape of the angular correlations, which is predominantly quadrupolar under general relativity. Here we explore a more flexible parameterization: we expand the angular correlations into a sum of Legendre polynomials and use a Bayesian analysis to constrain their coefficients with the 15 yr pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). When including Legendre polynomials with multipolesℓ≥ 2, we only find a significant signal in the quadrupole with an amplitude consistent with general relativity and nonzero at the ∼95% confidence level and a Bayes factor of 200. When we include multipolesℓ≤ 1, the Bayes factor evidence for quadrupole correlations decreases by more than an order of magnitude due to evidence for a monopolar signal at approximately 4 nHz, which has also been noted in previous analyses of the NANOGrav 15 yr data. Further work needs to be done in order to better characterize the properties of this monopolar signal and its effect on the evidence for quadrupolar angular correlations. 

Abstract Pulsar timing arrays (PTAs) are designed to detect low-frequency gravitational waves (GWs). GWs induce achromatic signals in PTA data, meaning that the timing delays do not depend on radio frequency. However, pulse arrival times are also affected by radio-frequency-dependent “chromatic” noise from sources such as dispersion measure (DM) and scattering delay variations. Furthermore, the characterization of GW signals may be influenced by the choice of chromatic noise model for each pulsar. To better understand this effect, we assess if and how different chromatic noise models affect the achromatic noise properties in each pulsar. The models we compare include existing DM models used by the North American Nanohertz Observatory for Gravitational waves (NANOGrav) and noise models used for the European PTA Data Release 2 (EPTA DR2). We perform this comparison using a subsample of six pulsars from the NANOGrav 15 yr data set, selecting the same six pulsars as from the EPTA DR2 six-pulsar data set. We find that the choice of chromatic noise model noticeably affects the achromatic noise properties of several pulsars. This is most dramatic for PSR J1713+0747, where the amplitude of its achromatic red noise lowers from ${\log }_{10}{A}_{\mathrm{RN}}=-{14.1}_{-0.1}^{+0.1}$to $-{14.7}_{-0.5}^{+0.3}$, and the spectral index broadens from ${\gamma }_{\mathrm{RN}}={2.6}_{-0.4}^{+0.5}$to ${\gamma }_{\mathrm{RN}}={3.5}_{-0.9}^{+1.2}$. We also compare each pulsar's noise properties with those inferred from the EPTA DR2, using the same models. From the discrepancies, we identify potential areas where the noise models could be improved. These results highlight the potential for custom chromatic noise models to improve PTA sensitivity to GWs. 

Abstract The millisecond pulsar J1713+0747 underwent a sudden and significant pulse shape change between 2021 April 16 and 17 (MJDs 59320 and 59321). Subsequently, the pulse shape gradually recovered over the course of several months. We report the results of continued multifrequency radio observations of the pulsar made using the Canadian Hydrogen Intensity Mapping Experiment and the 100 m Green Bank Telescope in a 3 yr period encompassing the shape change event, between 2020 February and 2023 February. As of 2023 February, the pulse shape had returned to a state similar to that seen before the event, but with measurable changes remaining. The amplitude of the shape change and the accompanying time-of-arrival residuals display a strong nonmonotonic dependence on radio frequency, demonstrating that the event is neither a glitch (the effects of which should be independent of radio frequency,ν) nor a change in dispersion measure alone (which would produce a delay proportional toν⁻²). However, it does bear some resemblance to the two previous “chromatic timing events” observed in J1713+0747, as well as to a similar event observed in PSR J1643−1224 in 2015. 

Abstract The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays (PTAs) through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$power law of the GWB strain spectrum. To do this, we create a semianalytic SMBHB population model, fit to North American Nanohertz Observatory for Gravitational Waves (NANOGrav’s) 15 yr GWB amplitude, and with 1000 realizations, we study the populations’ characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power law. The first, at 2 nHz, is below our GWB realizations with ap-value significancep= 0.05–0.06 (≈1.8σ–1.9σ). The second, at 16 nHz, is above our GWB realizations withp= 0.04–0.15 (≈1.4σ–2.1σ). We explore the properties of a loud SMBHB that could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by 3 orders of magnitude, from ∼10⁶to ∼10³, between 2 and 20 nHz. This causes a break in the strain spectrum as the stochasticity of the background breaks down at ${26}_{-19}^{+28}\mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the 26 nHz break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early Universe.